Muon-spin relaxation measurements on the dimerized 
spin-1/2 chains NaTiSi 2 6 and TiOCl 
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We report muon spin relaxation (/iSR) and magnetic susceptibility investigations of two Ti 3+ 
chain compounds which each exhibit a spin gap at low temperature, NaTiSi2C>6 and TiOCl. From 
these we conclude that the spin gap in NaTiSi2 06 is temperature independent, with a value of 
2 A = 660(50) K, arising from orbital ordering at Too = 210 K; the associated structural fluctuations 
activate the muon spin relaxation rate up to temperatures above 270 K. In TiOCl we find thermally 
activated spin fluctuations corresponding to a spin gap 2A = 420(40) K below T c \ = 67 K. We also 
compare the methods used to extract the spin gap and the concentration of free spins within the 
samples from ^iSR and magnetic susceptibility data. 

PACS numbers: 76.75.+i, 75.50.Ee, 75.10.Jm, 75.40.Cx 



The interplay between spin, charge, and orbital degrees 
of freedom is particularly subtle in strongly-correlated ox- 
ides containing octahedrally-coordinated Ti 3+ ions {t\ g ). 
Such compounds are typically insulators because the 
5 = 1/2 spins are localized in the t\ g orbitals. In ox- 
ides containing chains of Ti 3+ ions, superexchange via 
oxygen gives rise to antifcrromagnetic coupling. A well 
known instability that can affect half-integer spin chains 
is the spin-Peierls (SP) transition 1 ", in which magnetoe- 
lastic coupling dimerizes the chain and allows a spin-gap, 
2 A, to open below a characteristic temperature, Tsp, re- 
sulting in a spin-singlet ground state, (e.g. Tsp = 14 K 
for CuGeOiA T SP = 18 K for MEM(TCNQ) 2 3 ) In con- 
trast, for integer spins a Haldane gap will be present, 
precluding any low temperature instability^. 

Recently, two oxides, NaTiSi 2 6 (NTSO) and TiOCl, 
have been intensively studied because they undergo 
dimerization transitions at unusually high temperatures. 
NTSO has the pyroxene structure with chains of TiOe 
octahedra that are only weakly coupled to one another—. 
TiOCl has TiO bilayers within the ab plane, well sepa- 
rated by Cl _ ions. At low temperature dimerized chains 
of Ti 3+ ions form along the b axis, with their spins cou- 
pled by direct exchange 6 . The magnetic susceptibility, \, 
drops sharply at 210 K for NTSO 5 . and 67 K for TiOCli 
due to the opening of a spin gap. However, for both com- 
pounds, the ratio of the spin gap (which we will quote 
in Kelvin) to the dimerization temperature is larger than 
for the canonical SP case (for which 2A/Tsp = 3.53) 1 , 
demonstrating that the dimerization transitions are more 
complex than the canonical SP case, and pointing to the 
possible role of orbital physics. 

X-ray diffraction studies of NTSO show that the Ti 3+ 
chains dimerize below 210 K. 8 Phonon anomalies mea- 
sured using Raman scattering^ are consistent with the 



dimerization being driven by an orbital ordering at 
Too = 210 K, below which the system is condensed 
in one of two possible orbitally ordered spin-singlet 
states, breaking translational symmetry. This model 
is supported by a number of theoretical studiesiSiii* 1 ^, 
although a composite 5 = 1 Haldane chain ground 
state, with spin-triplet dimers of ferromagnetically cou- 
pled Ti 3+ spins, has also been proposed 13 . (See also 
Refs.[H0) 

Unusually, the dimerization of TiOCl occurs in two 
stages. There is a second order phase transition at 
T C 2 = 91 K into an incommensurate dimerized phase, and 
a first order transition at T c \ = 67 K into a commensu- 
rate SP phase&LiSiiSiii. Although early work suggested 
that orbital fluctuations may play a dominant role in 
driving the transitioni&i&iSiSI, recent experiments have 
ruled this outiiSi* 2 ^. In particular, optical measure- 
ments, in combination with a cluster calculation, revealed 
that the crystal field splitting is large enough to quench 
the orbital degree of freedomil. It appears that frus- 
tration between the two staggered chains in the bilayer 
crystal structure of TiOCl leads to an incommensurate 
SP state between T cl and T e2 15 i 23 i 24 . The magnitude of 
the spin gap has previously been measured using opti- 
cal techniques 1 ^ 5 , and NMRi 5 .. The former suggest a 
value of 2A ~ 430 K, whereas the latter gave a value of 
A = 430 ± 60 K, which is very large compared to the 
observed transition temperatures. 

Muon spin relaxation (/xSR) experiments 26 probe 
magnetic ordering and dynamics from a microscopic 
viewpoint. The applicability of this technique to 
spin-gapped systems has been demonstrated by stud- 
ies of the canonical SP compounds CuGeOjSLS, and 
MEM(TCNQ) 2 29 ' 30 , and the Haldane chain compound 
Y2BaNiOg2i. In this paper we present the results of /iSR 



experiments on NTSO and TiOCl. We also made mag- 
netic susceptibility measurements using a SQUID mag- 
netometer that allow us to compare our results with pre- 
vious datck^iii. We have measured the magnitude of the 
spin gap and found the concentration of unpaired spins 
in the dimerized state. 

Our polycrystalline sample of NTSO was prepared by 
a solid-state reaction^ of Na 2 TiSi40n, Ti, and Ti02- 
The sample of TiOCl was composed of small single crys- 
tals synthesized using standard vapour-transport tech- 
niques^ from Ti02 and TiCl3. Our ^lSH experiments 
were carried out using the MuSR and ARGUS spec- 
trometers at the ISIS facility, United Kingdom. These 
were done in zero applied magnetic field (ZF) and in 
small magnetic fields along the axis of the initial muon 
spin polarization. Spin-polarized positive muons (/i + , 
mean lifetime 2.2 /j,s, momentum 28MeV/c, 7 M = 2ir x 
135.5 MHzT -1 ) were implanted into the polycrystalline 
samples where they stop within ~ 1 ns. The decay 
positron asymmetry function 26 , A(t), is proportional 
to the average spin polarization of the muons stopped 
within the sample. Examples of the measured asymme- 
try spectra are presented in Fig. ^ 

The absence of coherent muon precession, which would 
be indicative of a spontaneous magnetic field, together 
with the observation that both high and low tempera- 
ture spectra relax to the same value with negligible miss- 
ing asymmetry (Fig. [Q, excludes the presence of long- 
range magnetic order in either material. The form of the 
muon spin relaxation is dependent on the distribution 
and time-dependence of the local magnetic fields around 
the site where the muon is implanted. In spin-gapped 
materials we write A(t) as a product of relaxation func- 
tions2L2£i2£. We approximate a Gaussian Kubo-Toyabe 
function^ 6 , (which models the Gaussian distribution of 
randomly orientated nuclear spins) using exp[— (at) ], 
use exp(— Xt) to describe the depolarization due to fluc- 
tuating electronic spins, and add a constant background 
Abg to describe those muons landing outside the sample; 
giving a fitting function: 

A(t) = A(0) exp[-((rf) 2 ] exp(-Ai) + A BG . (1) 

At high temperature the relaxation due to the fluctua- 
tions of nuclear spins dominates. Because the observed 
value of a is small we do not see a recovery in the muon 
asymmetry at longer times; this is also due, in part, to the 
presence of the electronic fluctuations. Fitting the data 
with (j as a free parameter over the whole temperature 
range, 15 — 340 K, showed that in NTSO it was, within 
experimental error, temperature independent and it was 
subsequently fixed at a — 0.06 MHz. In TiOCl the situa- 
tion is more complicated^. The values of a obtained for 
TiOCl are shown in the inset to Fig.[3Ja). 

The relaxation of the muon spin due to the fluctuat- 
ing electronic spins changes significantly as a function of 
temperature, dominating the relaxation at low tempera- 
ture. Providing the electronic spin-fluctuation rate v is 

fast compared with w(-B 2 )/7^, where is the fluctu- 
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FIG. 1: (Color online) Examples of muon decay asymmetry 
data in (a) NaTiSi 2 6 and (b) TiOCl, with fits to Eq.0 



ating field at the muon-site due to the electronic spins 
(the fast-fluctuation regime), then^S A cx 1/v. The tem- 
perature variation of A is shown in Fig.|2a) (NTSO) and 
Fig. EJa) (TiOCl). At low temperature, one possible re- 
laxation mechanism for the muon spin is via thermally ac- 
tivated electronic spins fluctuating across the spin gap2^. 
This would lead to the relaxation rate A(T): 

A(T) = Acxp(2A/T), (2) 

where A is a constant. X-ray measurements show that 
the degree of dimerization is essentially constant below 
Too (NTSO) 8 and T cl (TiOCl)^, so we take A as a con- 
stant for each compound. If an additional temperature- 
independent relaxation mechanism is present then: 

A(T)=A /[1 + Sexp(-2A/T)], (3) 

where B is a constant. 

For NTSO, Eq. 0fits the data quite well at tempera- 
tures up to Too and even slightly above, suggesting that 
the relevant spin fluctuations in the paramagnetic state 
are controlled by a similar energy scale (in fact, the an- 
tiferromagnetic exchange J ~ A)ii. There is no sharp 
change in A at Too? as might be expected if Haldane 
chains 1 ™ formed at that temperature, since the local field 
at the muon site should change significantly in that case. 
The best fit to the data, found by including points up to 
270 K, gives 2A = 660 ± 60 K. 

In contrast, for TiOCl, the behavior is more complex. 
We see no significant changes in A either at T C 2 or at 
T* ~ 120 K, the temperature at which a pseudogap has 
been proposedi^ii&S^. However, A rises sharply below 
T c i and we are able to fit the A values between 58 K 
and 68.2 K to Eq. H yielding 2A = 420 ± 40 K. This 
value is in agreement with the magnitude of the spin 
gap previously measured using other techniquesi^i^i 2 ^^. 
At lower temperatures, A shows a much slower increase 
with decreasing temperature, and Eq. ^only poorly de- 
scribes the data. Instead, a square root exponential, 
exp(— A/At), as would be expected for a dilute distribu- 
tion of slowly fluctuating electronic moments^ 6 -, provides 
a better description of the data. This interpretation is 
supported by longitudinal-field measurements which im- 
ply that the fluctuations leave the fast fluctuation regime 
when A > 0.16 MHz (In contrast, similar measurements 
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FIG. 2: (Color online) Fitted data for NaTiSi 2 6 . (a) Elec- 
tronic relaxation rate, A (Eq.0. The red solid line is a fit to 
Eq. |21 (b) Magnetic susceptibility, Xm , after subtraction of 
the low-temperature Curie Weiss tail, with the temperature- 
independent paramagnetic response marked xo- The blue line 
is a fit to Eq. El 

suggest that NTSO remains within the fast fluctuation 
regime at all measured temperatures). At much lower 
temperatures, we speculate that the dilute distribution 
of electronic moments originates from unpaired spins due 
to defects and/or impurities. 

To gain a rough estimate of n, the concentration of 
unpaired spins within the samples, from the /xSR data, 
we followed the method used in Rcf. 3(1 We assume the 
unpaired spins to have S = 1/2, given the spin of the 
Ti 3+ ions, although a defect in a chain would produce 
two free spins and any impurity creating the defect may 
itself have a spin, so our value will be larger than the con- 
centration of defects and impurities. Measurements were 
made in small magnetic fields, < 5mT for NTSO and 
2 mT and 10 mT for TiOCl, applied along the direction 
of the initial muon spin polarization, at low temperature. 
The asymmetry spectra were fitted using the product of 
a longitudinal field Kubo-Toyabe function^ and an ex- 
ponential relaxation to model the weak dynamics^. Fit- 
ting the data to this function we can extract the width 
of the distribution of local fields at the muon site. From 
this we estimate the concentration of unpaired spins sur- 
rounding the muon using the expression given in Ref . |37] 
and adapted to the muon case in Ref. |3(J For NTSO 
this gives an impurity concentration n = 1.7(3)%, and 
for TiOCl we find n = 1.1(2) %. 

Magnetic susceptibility, x, measurements (Fig. IHb) 
and Fig. G^b)) showed that our samples were compa- 
rable with those used in previous studies. From the 
low temperature Curie tail in x we estimate the con- 
centration of unpaired spins to be n = 2.10(4) % for 
NTSO and 0.6(1) % for TiOCl. The close agreement 
between these values and those obtained from /iSR sup- 



FIG. 3: (Color online) Fitted data for TiOCl. (a) Electronic 
relaxation rate, A (Eq. 0. The red solid line is a fit to 
Eq.|5|between 58K and 68. 2K, where the relaxation rate was 
found to be thermally activated. (Inset) Nuclear relaxation 
rate, a, (Eq. showing a weak temperature dependence' 3 -, 
(b) Magnetic susceptibility, Xm, after subtraction of the low- 
temperature Curie-Weiss tail. The solid blue line is the fit, 
described in the text, using the values of the dimerization pa- 
rameter from Ref. [23 to define the temperature variation of 
the spin gap and applying Eq.|l]to describe the susceptibility. 

ports our model describing the muon depolarization. The 
temperature-independent contribution to the susceptibil- 
ity was xo = 1-21(1) x 10- 9 m 3 mor 1 for NTSO, consis- 
tent with previous data-. For TiOCl we found that xo 
was negligible, once the contribution to the susceptibil- 
ity from the dimerization was considered. To extract 
the size of the spin gap from the magnetic susceptibil- 
ity data, the Curie and temperature-independent terms 
were subtracted and we take the contribution, Xm, due 
to the thermal activation of spins across the spin gap of 
magnitude 2A to be: 

X m = Cexp(-2A/fc B T). (4) 

The result of fitting Eq.0]to the data for NTSO is shown 
in Fig. [21b). Using a constant value for the spin gap 
in NTSO, this gives an excellent parameterization of 
the data from low temperature to 200 K and leads to 
a value of the spin gap of 2A = 595(7) K, somewhat 
smaller than the /^SR measurements suggest. The situa- 
tion in TiOCl is more complicated, since the dimerization 
varies strongly with temperature. The structural dimer- 
ization gives alternating exchange constants along the 
chain J1.2 = J(l ± 5) and we take S to vary linearly with 
the structural dimerization—. To model the form of the 
susceptibility we used the intensity of the superstructure 
reflection observed in X-ray diffraction measurements^ 3 *. 
For small values of 6 the intensity of this peak varies as 
S 2 . The spin gap is proportional 3 ^ to 5 2 ^ 3 and, using 
these results to convert the peak intensities into an ef- 
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TABLE I: Parameters derived from /xSR and magnetic sus- 
ceptibility data, x- The s Pi n g a P> 2A, the concentration of un- 
paired spins, n, and T m f = 2A/3.53, the temperature at which 
mean-field theory estimates that the SP transition should oc- 
cur, given the value of 2A measured by jj.SK. 



Sample NTSO TiOCl 

2A (K) /iSR 660(60) 420(40) 

2A (K) x 595(7) 
n (%) mSR 1.7(3) 1.1(2) 

n (%) x 2.10(4) 0.6(1) 

T mf (K) 190(20) 120(10) 



fective spin gap, we then used Eq. 01 to derive the line 
plotted in Fig. Etb). This was scaled by the spin gap 
value derived from the /zSR measurements and the con- 
stant C was chosen to fit the data, given we have no a 
priori values to choose for either case. A Curie- Weiss 
term was added to give the best fit to the low tempera- 
ture data, a process complicated by oxygen adsorption, 
and so there remains a small discrepancy between the 
model and the data. Given that this model should only 
work in the limit A 3> it can be seen to be remark- 
ably successful in describing the form of the data below 
80 K, showing that the susceptibility is indeed varying 
with the structural dimerization. The breakdown of this 
model may explain why magnetic measurements suggest 
a lower value of T c i than the XRD measurements. 
In conclusion, we have measured the magnitude of the 



spin gap using /iSR in the two Ti 3+ chain compounds 
NaTiSi20g and TiOCl. Both the magnetic susceptibility 
and /iSR data for NaTiSi2 0g are well described by as- 
suming that the spin gap is independent of temperature 
below the transition at Too = 210 K, and we find the 
magnitude of the spin gap, 2A = 660 ± 60 K, to be con- 
sistent with magnetic susceptibility measurements (Ta- 
ble [l| and recent theoretical workSiifliii. TiOCl shows 
a thermally activated increase in the muon relaxation 
rate below T cl = 67 K. Our value of 2A = 420 ± 40 K 
is consistent with previous NMR data— ^ and optical 
measurementsi&S^. Using mean-field theory to compare 
the measured spin gap in TiOCl to the observed transi- 
tions, we note (Table that T m f ~ T*, the temperature 
where psuedogap formation has been observecU^iS^, 
suggesting that T* is analogous to T sp in TiOCl, yet the 
system dimerizes at a far lower temperature. It is likely 
that the frustration between chains in the bilayer struc- 
ture lowers the temperature at which dimerization oc- 
curs. In TiOCl we can directly relate previously reported 
superstructure reflections^ to the magnetic susceptibil- 
ity. Comparing the concentrations of unpaired spins de- 
termined by /xSR and magnetic susceptibility measure- 
ments (Table we find that for both compounds the 
two methods are consistent, although our assumptions 
suggest that both techniques overestimate the impurity 
concentration. 

We are grateful to P. J. C. King and the staff of the 
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